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ABSTRACT 

A turbulent boundary layer flow analysis of the 
momentum and thermal flow fields near the forward stag- 
nation point due to a circular cylinder in turbulent 
crossflow is presented. Turbulence modeling length 
scale, anisotropic turbulence initial profiles and 
boundary conditions were identified as functions of 
the crossflow turbulence intensity and the boundary 
layer flow far field velocity. These parameters were 
used in a numerical computational procedure to calcu- 
late the mean velocity, .mean temperature, and turbu- 
lence double correlation profiles within the flow 
field. The effects of the crossflow turbulence on the 
stagnation region momentum and thermal flow fields 
were investigated. This analysis predicted the exist- 
ing measurements of the stagnation region mean velocity 
and surface heat transfer rate with crossflow Reynolds 
number and turbulence intensity less than 250 000 and 
0.05, respectively. 

NOMENCLATURE 

A 1 coefficients of Eq. (22), i = 1 to 7 
a- turbulence amplification factor, i = 1, 2, and 3 
B.j coefficients of Eq. (23), l = 1 to 4 

c correlation factor for the Reynolds stresses 

D cylinder diameter 

E sum of the dimensionless Reynolds normal stresses 

E i dimensionless Reynolds normal stress, i = 1, 2, 
and 3 

F dimensionless mean flow streamwise velocity 

f.j parametrical function of *, i = 1, 2, and 3 
G flow variable 

H dimensionless mean temperature 

i turbulence modeling length scale 

Nu Nusselt number 


Pr Prandtl number 

I 

Re d Reynolds number, D/v^ 

S dimensionless Reynolds shear stress, -uv/L^ 

T mean temperature i 

Tu crossflow turbulence intensity.y v^// 

t temperature fluctuation 

U mean flow velocity component along the X^ 

direction 

u velocity fluctuation along the Xj direction 

V transformed velocity of V 

V' mean flow velocity component along the X,, 

direction 

v velocity fluctuation along the ^ direction 

w velocity fluctuation along the Xj direction 

X 1 coordinates in the physical plane, l = 1 , 2 

and 3 

S stagnation point laminar boundary layer thick- 

ness 

n transformed coordinate along X^ 

e dimensionless variable, — "vt7 U 0 T e 

v kinematic viscosity 

t transformed coordinate along Xj 

? angle from the stagnation point 

<p Q the initial station 

— time mean average 

Subscripts 

e far field boundary condition 

oe far field boundary condition at the initial 

station 

1 m M grid point in the n direction 


1 


p pivot point or previous t station 

w surface condition 

» cross flow condition 


INTRODUCTION 

The momentum and thermal flow fields near the 
forward stagnation point of a circular cylinder in 
turbulent crossflow have been the focus of considerable 
study, (1_ to 1^) . One objective has been to understand 
the enchancement of the surface heat transfer due to 
crossflow turbulence near the leading edge of a fluid 
system component. Existing experiments, (3 and 4J, 
have shown that the local Nusselt number increased by 
75 percent as compared with laminar flow conditions 
when the turbulence intensity was 3 to 5 percent. In 
the initial turbine design stage, the blade leading 
edge is often modeled as a circular cylinder in a tur- 
bulent crossflow. With the current high operational 
gas temperature, the turbulence in the gas stream sig- 
nificantly increases the leading edge heat transfer. 
Thus, prediction of the stagnation region surface heat 
transfer becomes vital to the turbine blade cooling 
design. To estimate this surface heat transfer rate, 
correlations derived from experimental data may be 
used. The alternative is to rely on theoretical anal- 
ysis and numerical computation to obtain the surface 
heat transfer rate. 

Many two-dimensional unsteady laminar boundary 
layer flow analyses have examined the flow field in 
the vicinity of a stagnation point with the main stream 
velocity fluctuating about a steady mean. This type of 
analysis predicted little or no effects of free stream 
oscillation on the surface heat transfer rate for a 
flat plate or circular cylinder. The experiments on 
flat plates agreed in general with these theoretical 
results. However, the experiments, (JL to8), indi- 
cated a large increase in the heat transfer rate near 
the forward stagnation point of a circular cylinder. 
Theoretical studies, (9 and 10), showed that vortex 
stretching caused a substantial increase in the surface 
he-at transfer rate on cylinders. Therefore, the phe- 
nomena associated with a circular cylinder in & turbu- 
lent flow appear to be much different from a flat plate 
boundary layer flow. As the turbulent stream moves 
towards the stagnation point, amplification (or atten- 
uation) of the Reynolds stresses occurs due to the flow 
deceleration. Turbulence amplification because of vor- 
tex stretching also occurs as the flow turns and moves 
away from the stagnation point. These phenomena in- 
fluence the thin boundary layer and, therefore, the 
surface heat transfer rate. 

Smith, (5), used two-dimensional boundary layer 
theory to analyze the flow at the stagnation point on a 
circular cylinder. The eddy viscosity was assumed to 
be proportional to the turbulence in the free stream. 
The proportionality constant was determined from the 
experimental data. Reynolds analogy, with a tuijbulent 
Prandtl number of unity, was also employed. This anal- 
ysis indicated a greater effect of turbulence on heat 
transfer than on skin friction. Similar analytical 
approach was used by Gorla, ( UJ , to investigate the 
effects of unsteady free stream and free stream turbu- 
lence on skin friction at a stagnation point. The free 
stream turbulence effect on skin friction was found to 
be less than that was reported in (5). These existing 
analyses did not study the anisotropic turbulence flow 
field. Hijikata, (12), performed a theoretical and 
experimental study of the stagnation point anisotropic 
turbulence. An equation for the anisotropic turbulence 
was added to the k-e turbulence model to analyze the 


measurements. Maximum Reynolds normal stress was found 
at a location far away from the surface' and” the normal 
stress decreased continuously in the direction normal 
to the surface. His analysis also showed that the 
turbulence length scale effect on the surface heat 
transfer rate was small compared with the effect of 
turbulence intensity. Donaldson, (1_3 and 14_), has 
developed a turbulent boundary layer flow 'analysis, the 
method of invariant modeling, to study boundary layer 
flow transition and generation of atmospheric-clear air 
turbulence. Using a turbulence modeling technique to 
form a closed system of conservation equations for tur- 
bulent flow field, the analysis can be used to investi- 
gate dynamics among the anisotropic turbulence and mean 
flow streamwise property gradients. Since large mean 
flow streamwise velocity gradient and anisotropic tur- 
bulence exist within the turbulent stagnation point 
flow field, the Reynolds analogy or the k-c turbu- 
lence model is not sufficient to analyze the flow 
field. Therefore, based on the method of invariant 
modeling, the present paper proposes a turbulent bound- 
ary layer flow analysis and a computational scheme to 
study the anisotropic turbulence flow field and its 
effect on the stagnation region surface heat transfer 
on a cylinder in turbulent flow. 

Reported herein are the theoretical analysis and 
numerical computational procedure used in this study. 

A summary of the conservation equations is presented. 
The turbulence modeling length scale is defined as a 
function of the turbulence. The initial profiles and 
far field boundary conditions are related to the cross- 
flow turbulence and the inviscid flow velocity. The 
conservation equations are written in finite difference 
forms. A computational scheme is developed and per- 
formed to simulate the existing experiments. The com- 
putational results are then discussed and compared with 
the existing measurements. 

THEORETICAL ANALYSIS OF THE FLOWFIELD 

The flow field of interest in the present study is 
a small region within 20° of stagnation point. A sche- 
matic of this flow region, together with the physical 
and transformed coordinate systems, are shown in 
Fig. 1. It is assumed that the mean flow is steady, 
incompressible two-dimensional turbulent boundary layer 
flow. The mean flow field is described by two- 
dimensional boundary layer flow continuity, momentum, 
and enthalpy equations. The momentum and enthalpy 
equations contain the turbulence double correlations, 
u7 and vT. These double correlations are related to 
other turbulence correlations and the two-dimensional 
mean flow field through the turbulence dynamic equa- 
tions. An analysis using the turbulence dynamic 
equations to keep track of the turbulence double cor- 
relations was developed in (_13 and _14.) • This type of 
turbulence analysis was employed in the present study. 
Following the details in (13 and 14), conservation 
equations to describe the TurbuleiiT stagnation point 
flow field were derived. This set of conservation 
equations requires turbulence modeling techniques to 
form a closed system of equations for the turbulent 
flow field investigation. Conservation equations and 
turbulence closure assumptions are given in (Jl_5) . The 
turbulence modeling length scale, coordinate transfor- 
mation, turbulent flow initial profiles and boundary 
conditions are pertinent to the results of the present 
investigation. They are described in the following 
sections. 

Turbulence Modeling Length Scale 

For the turbulent flat plate boundary layer flow 
analysis, the turbulent mixing length is usually 
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assumed to vary linearly in a small distance normal to 
the surface and then to remain constant, approximately 
0.1 of the boundary layer thickness, across the rest 
of the boundary layer. In the present study, turbu- 
lence modeling length scale, with similar profile in - 
the X? direction as the mixing length in a flat 
plate boundary layer, was used to model the higher 
order correlation terms (tendency-towards-isotropy 
terms, turbulence diffusion terms, and turbulence dis- 
sipation terms) in the turbulence dynamic equations. 

It was also assumed that increasing the crossflow tur- 
bulence intensity would increase the length scale near 
the surface. Thus, the pivot point location Xg p , 
moves towards the surface as crossflow turbulencl in- 
creases. In the outer region, the length scale, s. e , 
is approximately 0.1 of the theoretical stagnation 
point laminar boundary layer thickness, a. This length 
scale profile is shown in Fig. 2. The proportional 
factor between X 2 p and the turbulence, Tu, was found 
by matching the present computational results with the 
existing measurements ( 1 _ to 8 ). 


Coordinate Transformation 

The coordinate transformation is to provide scale 
factors in the Xi and X 2 directions so that, in the 
numerical computation, the Xj step size can be in- 
creased and the number of steps in the Xg direction 
can be reduced. Using the transformation, a constant 
step size in X 2 can be used and the rate of change of 
the dependent variables in the X^ direction is also 
much reduced. 

The following similarity coordinate transforma- 
tions were used here: 


Energy equation: 
2^F aH/ at + VaH/3n 


= - 2c y- aTg/ac 4 p-p a 2 H/an 2 + ( 2 1 ) a e / a n 
e 

Turbulence equations: 

2cFaEj/ac + VaE^/an 

4 EjF U g q .5 

= - 2c — |j — aUg/ac 4 — a (tE " aEj/an)/an 


4 a 2 E 1 / an 2 + 2S(2c ) 0,5 aF/a n 


2 c 


vE 


0.5 (E - 3E-, ) 


4 2cE, (- 2aF/ a n - rr aF/ a n • an /aX, - 

1 u e Ugi 


2 cF aE 2 / at + VaE 2 /3n 
= 5^-3 (iE 0,S aE 2 /a n )/an + a 2 E 2 /a n 2 


( 3 ) 


(4) 


5 (x x ) = f 1 (Ug/vJdXj and nfXj.Xg) = UgX 2 /v/(2c ) 0 ' 5 

An ( external far field mean flow velocity of 
U e = 2Vl sin tp was assumed. The dependent 
variables were also nondimensionalized: 


vE 0.5 (E - 3E ? ) 

‘T? 3 

2 

4 2cE ? (2aF/ac 4 p- aF/a n • an/aX, - (5) 

6 u e U„t 


F = U/Ug, H = T/T e , S = - uv/U 2 , e = - vt/UgTg, 


E : = ^/U 2 , E 2 = 7/U 2 , and £ 3 . 7/U 2 . 

Therefore, E = Ei 4 E 2 + E 3 . 

Through mathematical manipulations, the conserva- 
tion equations, (15), can be written in the following 
forms in the transformed coordinate system: 

Continuity equation: 


2 tF aE 3 / a c 4 VaEj/an 


2E 3 F 


U„ 


- 2c — g — aUg/ac 4 a (zE aE 3 /an)/an 


2 2 vE 0 - 5 (E - 3 E 3 } 

4 a 2 E,/3n 2 4 2c ^ 1 — — 


2 rF 

^ ct 3 T7 

Ug^ 


2caF/ac 4 aV/an 4 F = 0 


(!) 2cFaS/ac 4 VaS/ 3 n 


where. 


U 2 t vF , ,,V 4 . 0.5 v- 

V = — jj — 3 n/aX, + (2c) fp 

e e 


Momentum equation: 
2 c F 3F / 3C 4 V3F/3n 


= 2 ^ 1 u~ F } all / ac 4 a 2 F/an 2 4 (2c) 0,5 aS/an 


(2) 


= - 2C aUg/ac 4 3 ^ 3 ( «.E° ' 5 aS/ an) / an 
e 

4 a 2 S/3n 2 4 (2C) 0 ' 5 E 2 3F/3n 
r0.5 , 2 

+ 2tS (_ Vr ' t - ? 5 

e U e i 


and, 


( 6 ) 


(7) 
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2c F 36 / at + Vse/3n 

= - 2c y 1 aT /sc + 3 ^ 3 (iE°- 5 3 e/ 3 n )/ 3 n 
e 

+ 3 2 e/3n 2 + (2C)° #5 E 2 sH/3n 

+ 2^e /s F / 3 C + -g— 3 n / 3 X ^ • 3 F / 3 n - — g — — ~T~%} ^ ® ) 

\ u e u e 4 I) r' 


3 H/ 3 c = 0 were imposed along the X 2 direction at 9o 
station. " ----- — 

4. The boundary conditions, F = 0, H = T w /Tp. 

V = 0, E 3 = 0, Eg = 0, E 3 = 0, S = 0, and e = 0, 
were also imposed at the surface (X 2 = 0). 

5. To determine the values of Ei >oe , Eg oe , 
^3,oe* an d S oe at Xg ss », 

U oe =2V : sin ’o < 14 > 

and 


Equations (4), (5), ( 6 ), and (7) show the rela- 
tionships among the Reynolds normal stresses, shear 
stress, and the mean flow velocity. Equations (3) and 
( 8 ) relate the turbulence effect on the temperature 
field and, thus, the surface heat transfer rate. With 
the assumptions of the length scale for the turbulence 
modeling, Eqs. (1) to ( 8 ) were solved numerically for 
the momentum and temperature fields. The initial flow 
profiles and boundary conditions are analyzed in the 
following sections. 

Initial Profiles 

Because the crossflow turbulence affects the flow 
field, the initial profiles can not be assumed. They 
were determined analytically from the governing 
Eqs. (1) to ( 8 ). Assumptions for the mean flow stream- 
wise gradients were made and the conservation equations 
were reduced to a set of ordinary differential equa- 
tions. These equations were solved for the property 
profiles at the initial location, 9 = <? 0 . The anal- 
ysis, with assumptions, is described here: 

1. At 9 = cp 0 and Xg r: «*, (n = 8 .), aF/Sn = ,0, 

3 H / 3n = 0 , 3 E 3 / 3 h = 0 , aEg/sn = 0 , 3E 3 /3n = 0, 3S/3n = 0, 
and se/an = 0 were assumed. In addition, F = 1, 

Hal, and e = 0 were imposed as boundary conditions 
at Xg « <=. 

2. To determine the Reynolds normal and shear 
stresses, Ei(<p,Xg) = EiUn, Xg) [Up(® 0 )/U e ( 9 )] 2 , and 
Eg( cp, Xg ) = E 2 U 0 , Xg) [U e ( <p p )/U e (<p) were assumed in 
the vicinity of the stagnation point. Therefore: 


3E L / 3C = E : (c 0 , x 2 ) ( 


and 


Re d % 


-2 


-) 


2E 2 /3C - EgU 0 , x g) — 7 - 

{Re d V 

were found for small ® 0 . Similarly, 

aE 3 /3c = e 3 (c 0 ,x 2 ) f x (® 0 ) 
3S/3C = s(c 0 ,x 2 ) f 2 <® 0 ) 
and 

ae/sc = e (C 0 ,X 2 ) f 3 (<p 0 ) 


(9) 

( 10 ) 

(ID 

( 12 ) 

(13) 


were assumed. 

3. The mean flow is symmetric about 9 = 0 plane. 
The mean flow streamwise velocity and temperature gra- 
dients are zero at 9 = 0 plane. These mean flow prop- 
erties were also assumed at locations near the stagna- 
tion point, 9 « 0. Therefore, aF/ac = 0 and 


(15) 


were assumed at Xg « ■». 

Applying the assumptions, (Eqs. (9) to (13)), to 
the governing Eqs. (4) to ( 8 ), the following results 
were found at 9 = ® 0 and Xg w » ( n = 8 .): 


t l,oe 


vE 1 - 5 

oe 

3Ut 
oe e 


vE 


Re d*o 


7 U 77 


oe 


( du e /dc) 0 - u 


357T 

oe 


(16) 


1 

oe e 


2v 

2 2 
U oe l e 


vE 


1.5 

oe 


t 2 ,oe 


3U z 
oe e 


-2 


V E 


Re 


+ U 


UTS - 

oe 


(17) 


d’o 


1 

oe e 


O ^ 

“2 — ? 
U 1 
oe e 


f l<< 


oe 


(du /de) 0 


vE 


TJ 


1.5 
oe 
I 


oe 3,oe 


r0.5 7 

1 uE oe 2 v 2 

3t U t ..2 2 
e oe e U Qe * e 


Eg(® 0 ) 

and. 


jj— (dU e /dc ) 0 - jp- 
oe oe 


f 3 (® 0 ) = 0 


p0.5 

t oe 


2 v 

..2 2 
U 1 
oe e 


(18) 


(19) 


( 20 ) 


where E oe = El oe + Eg oe + E 3 oe . Substituting 
Eqs. (14) and (15) into’Eqs. (16) to (19), Ej 0 p, Eg oe , 
E 3 oe , f 1 ( <p 0 ) > and fg(® 0 ) were determined. iJith tn£ 
assumptions in 2 , the above parameters were used to 
define the mean flow streamwise turbulence gradient 
terms except 3 S/ 3 C. 

To determine the 3 S/ 3 C at the initial location 
®o, the following relations among the Reynolds 
stresses, i.e.. 


S oe = c ( E l,oe E 2 ,oe^ 


1/2 


( 21 ) 
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at X2 ~ « was assumed. This relation is similar to 
that has been proposed for the flat plate boundary 
layer turbulence, ( 16J . For the present study, the 
factor, c = 0.0001, was used. 

Following the above detail, a set of ordinary 
differential equations with specified boundary condi- 
tions were established at <p 0 location. ; These equa- 
tions were solved numerically for the initial profiles 
at 9 0 = 0.0375° to avoid the singularity at the stag- 
nation point. 

Boundary Conditions 

Surface and far field boundary conditions are re- 
quired to analyze the downstream flow field. Phenome- 
nological boundary conditions for Ej e , E2ei and E3e 
were inferred from the results of existing heat trans- 
fer measurements. 

The surface heat transfer rate beneath a stagna- 
tion point flow field with given Rep and Tu is 
approximately constant in <?£20*. This surface heat 
transfer rate can be computed from the initial profile 

T ~T 

analysis. The Reynolds normal stresses, u and v 
and the cylinder axial direction Reynolds normal 

stress, \> , were derived from this heat transfer phe- 
nomenon. Numerical computation was performed to obtain 
satisfactory distributions of these Reynolds normal 
stresses to enable the numerical analysis to predict 
the surface heat transfer within ? <_ 20°. The numeri- 
cal method is described in the following Numerical 
Computation section. A rapid increase in the Reynolds 
normal stresses was necessary to obtain numerically the 
constant surface heat transfer rate near (p 0 . The 
rapid increase in the Reynolds normal stresses was con- 
sistent with the turbulence amplification along an 
accelerated streamline. 

Summary of the boundary conditions is given in 
the following: 

1. The turbulence a mplifica tion/att enuat ion fac- 
tors ar e de fined as a, = u 2 /uj|„, a 0 = . and 

p o 1 “ oe c e oe 

a, = vr/w £ These factors are assumed to be only 

functions of the mean flow streamwise velocity. A set 
of ai, a2 and 83 is presented in Fig. 3 together 
with the numerical results of surface heat transfer for 
Red = 1.85xl0 5 and T u = 0.03. 

2. The far field^ean flow streamwise velocity is 
described by U e a 2 vl sin 9. 

The far field Reynolds normal stresses are: 


E 


E 


E 


le 

2e 

3e 


7i u 2 

e e 


e e 


e e 


a. E. (U / U ) c 
1 l,oe oe e 


*2 E 2,oe <W‘ 


*3 e 3 .« <W 


while Ei 0 p, E2 oe and E3 oe were previously 
determined from {fie initial profile analysis. 

3. The following conditions, F e = 1, = 1, 

e=0, and S = c (Ei e • E2e) 0,B were also imposed at 
the far field boundary, p = 8. 

4. The surface boundary conditions at n = 0, are 
F = 0, V = 0, H = T w /T e , Ei = 0, E 2 = 0, E 3 = 0, S = 0, 
and 8=0. 

With the initial profiles at <p 0 and the above 
boundary conditions, the downstream flow field proper- 
ties were found by a downstream marching numerical 


computational procedure to solve the finite difference 
equations corresponding to the governing Eqs. (1) to 
( 8 ). 

NUMERICAL COMPUTATION 

The governing equations were written in finite 
difference forms to compute numerically the dependent 
variable distributions along the n direction at each 
t (or if) station. The following finite difference 
scheme and computational procedure were used. 

Finite Difference Scheme 

Every one of the governing Eqs. (1) to (8) can be 
expressed in a general form as: 

AjSG/h + A 2 8G/8n 

= A 3 G + A 4 3 ( Aj3G/ 8n) / an + A g a 2 G/an 2 + A ? (22) 

where A-j ( 1 = 1 to 7) may be constants, functions of 
the fluid properties and their gradients, or the trans- 
formation coordinates. 

For the initial profile calculation, the assump- 
tions and results described in the initial profile sec- 
tion were used to define aG/ ac - The finite difference 
approximation, aG/as = (G - G p )/at, where G p is the 
variable at the previous c station, was employed 
for downstream computation. The following finite dif- 
ference expressions were also used to approximate the 
derivatives along the n direction. 

(aG/an) m = (G m +} - G m _3)/(2*n) 

(a 2 G/ap 2 ) m * (Gm+i - 2G m + G m _i)/(an) 2 , 


and 

[a(A 5 aG/a n ) /an]m = [^^*112^+1 - ( A 5,m+l/2 + Vl/2) 
G m + ^.m-l/EGm-ll/Un)? 

with 


A 5,m+l/2 = ( A 5,m + 1 + A 5,m)/2 

A 5,m-l/2 = ( A 5,m + A 5,m-l)/2 

Substituting the above expressions into Eq. (22), the 
following form of the finite difference equation: 


B l G m + l + B 2 G m 


B 3 G m-l = B 4 


(23) 


was obtained. The coefficients B, (i = 1 to 4) are 
functions of A^ at node points m + 1, m, and m - 1. 
With the prescribed boundary conditions, the variable 
G at each node point m was calculated from Eq. (23) by 
simple successive-substitution method. 


Computational Procedure 

A trial-and-error iteration procedure was used. 

A flow chart for this numerical process was given in 
Fig. 4. The present computational procedures contain 
two main routines at each t station. First, the 
turbulent velocity flow field was calculated. This 
velocity flow field was then used to compute the cor- 
responding temperature distributions. Similar numeri- 
cal computation procedures were repeated at the next 
downstream location. 

At each 5 station, the coefficients A, were 
determined first. The variable profiles at the previ- 
ous t station was used in the first iteration and the 
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updated values at the current c station were then 
used to determine-the-coeff-icients for the next i-ter-a— - 
tion. However, the laminar stagnation point boundary 
layer profiles were used to compute the in the 
first iteration of the initial profile calculation. 

For the cases calculated in this^study, the compu- 
tational domain was limited to 0.0375° <_ 9 <_ 20° and 
n £ 8 . Constant grid dimension (an = 0.05) in n 
direction and variable stepsize along the c direc- 
tion (ac corresponding to a? = 0.0375° for 9 <_ 0.2° 
and a? = 0 . 2 ° for 9 > 0 . 2 °) were used. 

RESULTS AND DISCUSSION 

The above numerical computational procedure was 
applied to calculate the turbulence and the mean flow 
properties within 20 ° of the forward stagnation point 
of a circular cylinder in turbulent crossflow. The 
crossflow mean temperature was 333 K. The cylinder 
diameter was 10.15 cm and its surface temperature was 
constant at 278 K. The Prandtl number was 0.7. The 
magnitudes of the mean crossflow velocity and its tur- 
bulence were prescribed to have approximately the same 
cross flow Reynolds number and turbulence intensity as 
they were reported in the existing experiments. There- 
fore, the present analytical results could be compared 
with the measurements. The present comput atio nal re- 
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suits were derived with the assumption of v = v . 

These numerical results are described and discusseS 
here. 

1. The mean flow streamwise velocity profiles at 
the initial station (9 = 0.0375°) with Rep = 2.50xl0 5 
and Tu <_ 0.05 are shown in Fig. 5. These profiles 
were also compared with the theoretical laminar bound- 
ary layer solution, ( lf[) . Crossflow turbulence in- 
creased the near surface mean streamwise velocity and 
decreased the mean velocity at locations far away from 
the surface. Existing velocity profile measurements, 
at 9 = 0 ° with Tu = 0.05, (17), were also shown in 
Fig. 5. General agreement couTJ be found from the 
comparison. 

Examples of the computational results of the 
Reynolds stress profiles at the initial station are 
shown in Fig. 6 . The Reynolds normal stresses increase 
and then decrease along the radial direction towards 
the surface. Near the surface, high E 2 Reynolds 
normal stress was calculated. These profiles are dif- 
ferent from the Reynolds stress distributions within a 
flat plate turbulent boundary layer, (16). 

2. The Nusselt number at the initial location was 
computed from the mean temperature profile. The numer- 
ical results with different crossflow Reynolds number 
are compared with existing experimental results, (5), 
in Fig. 7. The correlation between the Nusselt num- 
ber, Reynolds number, and turbulence intensity from 
several existing experiments and the present analysis 
were also plotted in Fig. 8 . The present analysis 
satisfactorily predicts the stagnation point surface 
heat transfer rate over the ranges of Rej < 2.5X10- 1 
and Tu <_ 0.05. Although lower surface heaT transfer 
rate, as - compared with the results in (j>), was cal- 
culated, the present analysis predicted the stagnation 
point heat transfer rate correlations, (5 to 8 ), 
especially at small and large values of T u (Red) 1 ' z . 

The numerical computational results also show that the 
parameter, Nu(Re,j) _ ®‘ 5 , approaches a maximum value 
( = 1.6) at large TuJRep)’' 2 . This result is consis- 
tent with the stagnation point heat transfer assump- 
tion used in turbine blade heat transfer design, (18) . 

3. The mean flow streamwise velocity profiles at 
downstream locations are shown in Fig. 9 for case with 
Rea = 1.85xl(P and Tu = 0.03. Large variation in 
this mean velocity component occurred in 9 <_ 0.15 . 


Although turbulence amplification was imposed along the 
outer boundary of the flow field, rapid mean' flow “ 
acceleration reduced the turbulence double correla- 
tions with a factor of 1/ul. This reduction in the 
outer edge boundary conditions decreases the Reynolds 
stresses within the flow field. Thus, the mean flow 
streamwise velocity distribution tended toward the lam- 
inar boundary layer p-o f ile as flow moves downstream. 

Examples of the computational results of the 
Reynolds normal stress profiles at downstream locations 
are shown in Fig. 10. High Reynolds normal stresses, 
especially the Ej component, appear near the sur- 
face. Profiles with large Ej and small E 2 near the 
surface were calculated at downstream locations near 
9 = 20°. These results are similar to the Reynolds 
normal stress profiles in flat plate turbulent boundary 
layer. They are significantly different from the ini- 
tial profiles in Fig. 6 . 

The Reynolds shear stress distributions at three 
downstream locations are given in Fig. 11. Maximum 
shear stress occurs near the surface and the location 
shifts towards the surface as flow moves downstream. 
These results also show that shear stress profile re- 
stores to the flat plate turbulent boundary layer pro- 
file at locations about 9 = 20 °. 

4. Profiles of the dimensionless mean tempera- 
ture, H, and the turbulence correlation, e, with 

Rey = 1.85x10^ and Tu = 0.03 are shown in Fig. 12. 
Mean temperature prof'le variation only occurs near the 
stagnation point. Comparison with the laminar boundary 
layer temperature profile is also made. Crossflow tur- 
bulence increases the near surface mean temperature and 
decreases the mean temperature in the far field of the 
boundary layer. The a profiles show a large surface 
normal gradient at the initial location (9 ■ 0.0375°). 

As it can be seen from Eq. ( 8 ), large E 2 near the 
stagnation point induces this turbulence gradient. 

This large turbulence gradient then increases the mean 
flow temperature near the surface and, therefore, the 
surface heat transfer rate. Small values in e was 
computed at downstream locations. The effect of the 
0 variation on the mean temperature is mostly limited 
to a region away from the surface. 

5. In the downstream flow calculation, the turbu- 
lence amplification/attenuation factors were derived 
from the numerical computation to obtain an approxi- 
mately^constant Nusselt number distribution within 

9 <_ 20° region. With these boundary conditions, a 
negative Reynolds stress gradient was computed along 
most parts of the outer boundary of the computational 
domain. This is possible due to the existence of the 
mean flow velocity gradients. This negative gradient 
was reduced to approximately zero gradient at 9 = 20 ° 
locations where the Reynolds stress profiles are simi- 
lar to those of the flat plate turbulent boundary 
layer. These computational results implied that the 
assumption of zero Reynolds stress gradients as far 
field boundary conditions is valid at locations 
{9 ^ 20 °) away from the stagnation point. 

6 . The above results are subjected to the assump- 
tions and constraints in the multi-equation Reynolds 
stress modeling to account for the anisotropic turbu- 
lence. The length scale of the crossflow turbulence 
could induce different turbulence intensity near the 
stagnation point, (17), and its surface heat transfer 
rate, ( 8 ). The crossflow turbulence length scale 
effect was not considered in the present study. 

CONCLUSION 

Turbulence length scale, anisotropic turbulence 
initial and boundary conditions for a theoretical and 
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numerical boundary layer flow analysis of the stagna- 
tion region momentum and thermal flow fields have been 
defined as functions of the incoming crossflow turbu- 
lence and mean flow velocity. Methods to define the 
length scale, initial and boundary conditions and their 
applications to calculate the turbulence double corre- 
lations, mean velocity and temperature, and the surface 
heat transfer around the stagnation point of a circular 
cylinder in a turbulent crossflow were described. From 
the comparisons with the existing measurements, the 
present analysis was found capable of predicting accu- 
rately the local surface heat transfer rate with the 
crossflow Reynolds number less than 2.5x10® and turbu- 
lence intensity up to 5 percent. 
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Figure 2. - Turbulence modeling length scale. 
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Figure 3. - Anisotropic turbulence boundary conditions and 
surface heat transfer correlation. Re^, 1.85x10^; Tu, 0.03. 





Figure 4. - Flow chart of the numerical computational procedures. 
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Figure 6. - Examples of Reynolds stress profiles at initial station, 
Re d , 1. 85x10'’; (p Q , 0.0375°; Tu, 0.03. 







NUSSELT NUMBER, Nu 



.01 .02 .03 .04 .05 


TURBULENCE INTENSITY, Tu 


Figure 7. - Comparisons between analytical and experimental 
near stagnation point heat transfer. 
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Figure 8. - Stagnation point heat transfer correlations from the 
analyses and experiments. 
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Figure 9. - Examples of mean flow streamwise velocity profiles at 
different locations. Re^, 1.85x10-’; Tu, 0.03, 
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Figure 10 . - Example of downstream Reynolds normal stress profiles. 
Re d , 1. 85 X 10 5 ; Tu, 0. 03. 
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Figure 12. - Examples of the mean temperature and turbulence profiles. 
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